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Condensation phenomena arise through a collective behaviour of particles. They are ob¬ 


served in both classical and quantum systems, ranging from the formation of traffic jams in 
mass transport models to the macroscopic occupation of the energetic ground state in ultra¬ 
cold bosonic gases (Bose-Einstein condensation). Recently, it has been shown that a driven 
and dissipative system of bosons may form multiple condensates. Which states become the 
condensates has, however, remained elusive thus far. The dynamics of this condensation are 
described by coupled birth-death processes, which also occur in evolutionary game theory. 
Here, we apply concepts from evolutionary game theory to explain the formation of multiple 
condensates in such driven-dissipative bosonic systems. We show that vanishing of relative 
entropy production determines their selection. The condensation proceeds exponentially fast, 
but the system never comes to rest. Instead, the occupation numbers of condensates may os¬ 
cillate, as we demonstrate for a rock-paper-scissors game of condensates. 

Condensation phenomena occur in a broad range of contexts in both classical and quantum 
systems. Networks such as the World-Wide-Web or the citation network perpetually grow by the 
addition of nodes or links and they evolve by rewiring. Over time, a finite fraction of the links of 
a network may be attached to particular nodes. These nodes become hubs and thereby dominate 
the dynamics of the whole network; they become condensate nodeP®. Condensation also occurs 
in models for the jamming of traffid® and in related mass transport models in which particles hop 
between sites on a latticd^I^. A condensate forms when a finite fraction of all particles aggregates 
into a cluster that dominates the total particle flow. Bose-Einstein condensation, on the other 
hand, is a quintessentially quantum mechanical phenomenon. When an equilibrated, dilute gas 
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of bosonic particles is cooled to a temperature near absolute zero, a finite fraetion of bosons may 
eondense into the energetie ground statd^^i^. Long-range phase eoherenee builds up and quantum 
physies beeomes manifest on the maeroseopie seald ^^ l *"^ !. 

In both the elassieal and the quantum meehanieal eontext, eondensation oeeurs when one or 
multiple states beeome maeroseopieally oeeupied (they beeome eondensates), whereas the other 
states beeome depletecfi^H. However, the physieal origins of eondensation in the above examples 
differ from eaeh other. Why and how eondensation arises in a partieular system remains a topie of 
general interest and vivid researeh. 

Here, we study eondensation in two systems from different fields of researeh: ineoherently 
driven-dissipative systems of non-interaeting bosons and evolutionary games of eompeting agents. 
As we show below, the physieal prineiple of vanishing entropy production governs the formation 
of eondensates in both of these systems. The entities that eonstitute the respeetive system shall be 
called partieles. They may be quantum or elassieal particles (bosons or agents). The dynamies of 
these partieles eventually lead to eondensation into partieular states (quantum states or strategies). 
Before deseribing the above two systems, we now introduee the mathematieal framework of our 
study. 


On an abstraet level, we eonsider a system of S (non-degenerate) states Ei,i = 1,..., S', 
eaeh of whieh is oeeupied by iVj > 0 indistinguishable partieles, see Fig. ??a. The eonfiguration 
of the system at time t is fully eharaeterized by the oeeupation numbers N = (A^i, N 2 ,..., Ns). 
This eonfiguration ehanges eontinuously in time due to the transition of partieles between states. 
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The total number of particles in this coupled birth-death process is conserved (N = 

We are interested in the probability P(N, t) of finding the system in configuration N at time t. 
The temporal evolution of the probability distribution P(N, t) is governed by the classical master 


equatiorlfdh^ 








( 1 ) 


where e* G denotes the unit vector in direction i (equal to one at index i, otherwise zero). The 
rate for the transition of particles from state Ej to Ei depends linearly on the number of particles 
in the departure and in the arrival state: 


with rate constant > 0 and constants Sij > 0. 


( 2 ) 


Condensation in this framework is understood as the macroscopic occupation of one or mul¬ 
tiple state^l21i^ We consider a state Ei as a condensate when the long-time average of the number 
of particles in this state scales linearly with the system size ((iVj)t ~ 0{N) for large t). Hence, a 
condensate harbours a finite fraction of the total number of particles for large systems {N 3> 1). 
We refer to a state as depleted when its average occupation number scales less than linearly with 
the system size. Therefore, the fraction of particles in a depleted state vanishes in the limit of large 
systems. 

Depending on the values of the rate constants r^, numerical simulation of equations ([T]) with 
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rates @ reveals that all states, multiple states, or only one state beeome eondensates when the 
partiele density N/S is large enough to deteet eondensatiorpl. Thus far, various questions about 
eondensation have remained elusive for the eoupled birth-death proeess defined by equation ([T]): 
Whieh of the states beeome eondensates? How does this seleetion of eondensates proeeed? Is it 
possible to eonstruet systems that eondense into a speeifie set of eondensates? 

In the following, we answer these questions by illuminating the physieal prineiple that gov¬ 
erns the formation of multiple eondensates on the leading order timeseale. We show that the 
vanishing of relative entropy produetion determines the seleetion of eondensates (see equations @ 
and Q below). We elaborate how eondensate seleetion is determined by the rate eonstants r^. The 
eondensation proeeeds exponentially fast into a dynamie, metastable steady state within whieh the 
oeeupation numbers of eondensates may oseillate. By applying our general results to systems with 
many states, we show that the interplay between eritieal properties of sueh networks of stated 
and dynamieally stable network motif^ determines the seleetion of eondensates. The results 
of our analysis apply to any system whose dynamies are deseribed by the eoupled birth-death 
proeesses ([1]) with rates Q. Before proeeeding to the mathematieal and numerieal analysis of 
eondensation in these proeesses, we now give a brief overview of sueh systems. 

Results 

Non-interacting bosons in driven-dissipative systems. The elassieal master equation ([T]) has 
reeently been derived by Vorberg et al. in the study of bosonie systems that are dissipative and 
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driven by external souree^. For a system of non-interaeting bosons that is weakly eoupled to a 
reservoir and driven by an external time-periodie foree (a so-ealled Floquet systemj^SfSl, one ean 
eliminate the reservoir degrees of freedom (Born and Markov approximation and the density 
matrix of the system beeomes diagonal (see the Supplement of the work of Vorberg et al.E%. The 
effeetive dynamies of the bosons beeome ineoherent and are eaptured on a maeroseopie level in 
terms of the eoupled birth-death proeesses ([T]) with rates Fj^j = rij{Ni + l)Nj (that is all Sij = 1 
in the rates @). These non-equilibrium setups may not only lead the bosons into a single, but also 
into multiple eondensate^. 

For the ineoherently driven-dissipative systems deseribed above, the state Ei denotes a time- 
dependent Floquet stateP^^^. The total rate F i^j for the transition of a boson from state Ej to E^ 
depends linearly on the number of bosons in the departure state (Nj) and the arrival state (Nt -\- 
1). The latter faetor stems from the indistinguishability of bosons and refleets their tendeney to 
eongregate. Although, we refer to equation ([T]) as a elassieal master equation and eoherenee does 
not build up, the quantum statisties of bosons is still eneoded in the funetional form of The 
rate eonstant rij is determined by mieroseopie properties of the system and the reservoir. 

Condensation in the above setup is to be distinguished from Bose-Einstein eondensation. 
Typieally, studies on Bose-Einstein eondensation foeus on the existenee of long-range phase eo¬ 
herenee in thermal equilibriurrpHS]^ kinetie formatioiJi^IlEZHlll, and the fragmentation of a eo- 
herent eondensate into multiple eondensates (for example when the equilibrium ground state is 
degenerates 'l l In eontrast, the elassieal birth-death proeesses O with rates @ deseribe eonden- 
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sation in bosonic systems that are externally driven by a continuing supply of energy, dissipate into 
the environment, and exhibit decoherence. 

Equations of type ([T]) may also arise in atomic physics and quantum optics and are known as 
Pauli master equationd^^HS] describe how the population of S non-degenerate energy levels 
changes over time when a system harbours N indistinguishable, non-interacting bosonic atoms. 
Such changes may occur by interactions with a radiation field that induces transitions between 
energy levels. A theoretical description of these transitions in terms of a Pauli master equation is 
appropriate if coherence is negligible. As in the previous example, the system then approaches 
a state in which some of the energy levels are macroscopically occupied (condensates) whereas 
others are depleted. More generally, whenever a rate constant rij governs the transition of a single 
boson from a state Ej to Et, the rates @ with Sij = 1 for all i and j apply if N non-interacting 
bosons are brought into the systerrP^. 

Strategy selection in evolutionary game theory. The classical master equation ([T]) also occurs in 
evolutionary game theory (EGT). Historically, EGT was developed to study evolutionary processes 
that are driven by selection and mutatioi^UlZl and seeks to identify optimal strategies for compet¬ 
itive interactions. Eor example, EGT has been applied in the study of the prominent “rock-paper- 
scissors” game, which was proposed as a facilitator of species coexistence and has inspired both 
experimental and theoretical researchP®^. Eurthermore, the “prisoner’s dilemma” game serves as 
a paradigmatic model to explore the evolution and maintenance of cooperatiorl^lllll. The interplay 
between non-linear and stochastic effects underlies the dynamics of such evolutionary gamed^^®. 
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In EGT, one typically considers a system of N interacting agents (classical particles) who 
repeatedly play one fixed strategy Ei out of the S possible choices Ei, E 2 ,..., Es. In each suc¬ 
ceeding interaction, the defeated agent adopts the strategy of its opponent. Since Nj agents playing 
strategy Ej can potentially be defeated by one of the Ni agents playing strategy E^, the rate of 
change is = rijNiNj. If an agent who plays Ej can also spontaneously mutate into an agent 
who plays E^ (with rate ^ij = VijSij), one recovers the classical master equation ([T]) with rates 

Thus, there exists a correspondence between condensation in incoherently driven-dissipative 
bosonic systems and strategy selection in evolutionary game theory: the transition of bosons be¬ 
tween states can be interpreted in terms of the interaction and mutation of agents employing evo¬ 
lutionary strategies. In effect, the states in an incoherently driven-dissipative setup play an evolu¬ 
tionary game and the winning states form the condensates. 

After having introduced the above examples, we now proceed with the mathematical and 
numerical analysis of the classical master equation ([T]). We show that the dynamics of condensa¬ 
tion change on two distinct timescales. At the leading order timescale, the dynamics are described 
by a set of nonlinearly coupled, ordinary differential equations (see equation @ below), which 
determine the states that become condensates. We identify these states by applying concepts from 
EGT. After an exposition of the physical principles that underlie the condensation dynamics, im¬ 
plications of our general results for incoherently driven-dissipative systems are discussed. 

The antisymmetric Lotka-Volterra equation. The total number of particles needed for conden¬ 
sation phenomena to occur is large (N S> 1). In order to detect macroscopic occupancies, it is also 
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assumed that the particle density N/S is large. Therefore, one may approximate the classical mas¬ 
ter equation ([T]) by a Langevin equation for the state concentrations Xi{t) = Ni{t)/N (details of 
the derivation are provided in Supplementary Note 1). Originally proposed for Brownian particles 
suspended in a liquid, the Langevin equation decomposes the dynamics of a sample trajectory of 
the random process into two contributions: into a deterministic drift and into noise stemming from 
the discreteness of particle numbers (“demographic fluctuations”). Both the demographic fluc¬ 
tuations and the contribution to the deterministic drift that corresponds to mutations in the EGT 
setting are suppressed by a small prefactor 1/N. Therefore, these terms change the dynamics only 
slowly. The deterministic drift that corresponds to interactions between agents is, however, not 
suppressed. It thus governs the dynamics to leading order. 


Hence, we find that the leading order dynamics of the condensation process ([T])-@ are de¬ 
scribed by the differential equations: 


—Xi = Xi{Ax.)i. 


(3) 


The matrix A is antisymmetric and encodes the effective transition rates between states (a^j = 
Tij — rji). The constants Sij that occur in the definition of the rates Q do not change the leading 
order dynamics, but they become relevant on subleading order timescales. 


We refer to equation @ as the antisymmetric Lotka-Volterra equation (ALVE). It provides 
a description of pairwise interactions that preserve the total number of particles. Therefore, the 
ALVE finds a broad range of applications in diverse fields of research, in addition to the aforemen¬ 
tioned condensation of bosons far from equilibrium. It was first studied by Volterrapf^in the context 
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of predator-prey oscillations in population biologjEOllIll] plasma physics, the ALVE describes 
the spectra of plasma oscillations (Langmuir waves and in chemical kinetics it captures the 
dynamics of bimolecular autocatalytic reaction^f^UlMl. in EGT, the ALVE is known as the replica¬ 
tor equation of zero-sum games such as the rock-paper-scissors gameE^HHUl. Table [^summarizes 
all of the above analogies. 

Despite the simple structure of the ALVE, it exhibits a rich and complex behaviour. In 
the following, we show how the mathematical analysis of the ALVE explains condensation into 
multiple states (condensate selection). To this end, we extend an approach for the analysis of the 
ALVE that was introduced in the context of EGlI^^llin!. 

Production of relative entropy and condensate selection. Our analysis starts from a theorem 
in linear programming theor)^^. Given an antisymmetric matrix A, it is always possible to find a 
vector c that fulfils the following conditions: its entries are positive for indices in 
and zero for indices inJ = {1,..., Sj — I, whereas the entries of Ac are zero for indices in I 
and negative for indices in / (Eig. ??b). Although several vectors c with these properties may 
exist, the index set / is unique and, thus, determined by the antisymmetric matrix A. Einding 
such a “condensate vector” c is crucial for the understanding of condensate selection and of the 
condensation dynamics. The condensate vector has the following physical interpretations. 

All condensate vectors yield fixed points of the ALVE Because of the antisymmetry 
of matrix A, a linear stability analysis of these fixed points does not yield insight into the global 
dynamics (Supplementary Note 1). However, global stability properties can be inferred by showing 
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that the relative entropy of a eondensate veetor to the state eoneentrations, 


D{c\ |x) = ^Ci log (ci/xi ), (4) 

i£l 

is a Lyapunov function (note that we do not consider the relative entropy of the state concentrations 
to the condensate vector, but define the relative entropy vice versa). The relative entropy Q de¬ 
creases with time and is bounded from below (see Methods and Supplementary Fig. 1). Therefore, 
the dynamics relax to a subsystem in which relative entropy production is zero. The relaxation of 
relative entropy production is reminiscent of Prigogine’s study of open systems in non-equilibrium 
thermodynamics. Indeed, we find that the system, to cite Prigogine’s phrase, “settles down to the 
state of least dissipation’^^. 

This state of least dissipation is characterized even further by the condensate vector c. Con¬ 
sidering the definition of the relative entropy Q and its boundedness, it follows that every concen¬ 
tration Xi with i G I remains larger than a positive constant. On the other hand, states with indices 
in I become depleted for long times (see Methods). Therefore, we find that the condensate vector 
determines the selection of condensates. Positive entries of c correspond to states that become 
condensates, whereas zero entries of c correspond to states that become depleted. Both the set 
of condensates and the set of depleted states are unique (Fig. ??b) and independent of the initial 
conditions. Generically, the entries of the condensate vector are also unique upon normalization 
(its entries sum up to one) and yield the rate |(^c)j| at which a state Ei becomes depleted. The 
condensate selection occurs exponentially fast (see Fig. ?? and Methods). 

After relaxation, the dynamics of the system are restricted to the condensates. In other 
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words, the condensates form the attractor of the dynamics. However, the dynamics in this sub¬ 
system do not come to rest. The state of least dissipation is a dynamic state with a perpetually 
changing number of particles in the condensates: periodic, quasiperiodic, and non-periodic oscil¬ 
lations are observed (Fig. ??b and Supplementary Fig. 1). In the generic case, the entries of the 
condensate vector represent the temporal average of condensate concentrations according to the 
ALVE Q. After condensate selection, the dynamics of these active condensates take place on a 
high-dimensional, deformed sphere!^. 

An algebraic algorithm to find the condensates. Numerical integration of the ALVE @ is 
neither a feasible nor a reliable method for identifying condensates (Fig. ??, Supplementary Figs. 1 
and 2). Instead, we determine these states by numerically searching for a condensate vector c. To 
this end, we reformulate the above conditions on c in terms of two linear inequalities 

Ac < 0 and c — Ac > 0 . (5) 

We solve these inequalities with a linear programming algorithm that is both reliable and efficient. 
The time to find a condensate vector scales only polynomially with the number of states S (see 
Supplementary Fig. 3 and Methods for details). 

Condensation in large random networks of states. We used our combined analytical and nu¬ 
merical approach to study how the connectivity of a random network of states affects the selection 
of condensates under the dynamics of the ALVE Q. The connectivity specifies the percentage of 
states between which particle transitions occuit^^!^. After having generated a network with a given 
connectivity, the strength and direction of an allowed transition between states Ej and Ei were 


12 




determined by randomly sampling the eorresponding effective rate constant Uij = Vij — rji. 

Our results for condensation in large random networks of states are summarized in Fig. ??. 
When the connectivity of a network is zero, all of its states are isolated. Particles are not ex¬ 
changed between states and none of the states becomes depleted. For an increased connectivity, 
isolated pairs of states are sampled in a random network. One state in an isolated pair is always 
depleted and the average number of condensates decreases rapidly. Upon approaching a critical 
connectivity, cycles and trees of all orders become embedded in a random network. This critical 
connectivity scales inversely with the number of stated^. We observe that, under the dynamics 
of the ALVE the average number of condensates becomes minimal for a connectivity that 
also scales inversely with the number of states (see Fig. ??g). We attribute this minimum to the 
interplay between the criticality of random networks and condensate selection on connected com¬ 
ponents of the network^n^. Embedded directed cycles are a recurring motil^ in the remaining 
network of condensates after condensate selection. Above the critical connectivity, a single giant 
cluster is formed. On average, half the number of states in this giant cluster become condensates 
once the network is fully connected (C = 1) (Fig. ??aj— I— 1^. Thus, our analysis emphasizes the 
importance of critical properties of random networks for condensate selection. 

Design of active condensates. Our understanding of condensate selection can be used to design 
systems that condense into a particular network of states, a game of condensates. We exemplify 
this procedure by formulating conditions under which a system relaxes into a rock-paper-scissors 
(RPS) game of condensate^. Three particular states Ei,E 2 , and E 3 in a system become a RPS 
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cycle of condensates if, and only if, the following two eonditions are fulfilled (Fig. ??). First, the 
“RPS eondition” requires that the rate eonstants between the three states form a RPS network (for 
example, ri 2 > r 2 i,r 23 > ^ 32 , and r 3 i > ri 3 ). Seeond, the “attraetivity eondition” requires that 
the inflow of partieles into the RPS eyele from any other state is greater than the outflow to that 
state Ek (for all A; = 4,..., S'). The values of the rate eonstants between the states that beeome de¬ 
pleted are irrelevant. More eomplex games of eondensates ean be designed by formulating similar 
eonditions on the rate eonstants. These eonditions are formulated as inequalities that depend on 
Pfaffians of the antisymmetrie matrix A and its submatriees (see Methods j^. The flow of partieles 
between states in these systems eauses eondensate eoneentrations to oseillate (Fig. ??b). 

Discussion 

Our findings thus suggest intriguing dynamies of eondensates in systems whose temporal evo¬ 
lutions are eaptured by the elassieal master equation ([T]) with rates ([^; for example in driven- 
dissipative systems of non-interaeting bosons. Condensates observed on the leading order timeseale 
are metastable. For longer times, relaxation into a steady state oeeur^l^J^. When detailed balanee 
is broken in the system of eondensates, the net probability eurrent between at least two states does 
not vanish and a non-equilibrium steady state is approaehed^^H^. The simplest way of designing 
sueh eondensates is illustrated by the above RPS game. In this game, detailed balanee is bro¬ 
ken, for example, when the transition of partieles is unidireetional (with totally asymmetrie rate 
eonstants ri 2 > r 2 i = 0,r23 > r 32 = 0, and r 3 i > ri 3 = 0). For non-interaeting bosons in driven- 
dissipative systems, the eontinuing supply with energy through the external time-periodie driving 
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force (Floquet system) and the dissipation of energy into the environment may, therefore, prevent 
the system from reaehing equilibrium. How sueh systems may be realized in an experiment poses 
an interesting question for future researeh. 

The transition of partieles between eondensates in the here studied eoupled birth-death pro- 
eesses parallels the interaetion and mutation of winning agents in evolutionary game theory, re- 
fleeting an “evolutionary game of eondensates”. Our results suggest the possibility of ereating 
novel bosonie systems with an oseillating oeeupation of eondensates. Non-interaeting bosons in 
ineoherently driven-dissipative systems are promising eandidates. Sinee the antisymmetrie Lotka- 
Volterra equation also arises in population biology, ehemieal kineties, and plasma physios, all of 
our mathematioal results apply to these fields as well. 

Methods 

Asymptotics of the antisymmetric Lotka-Volterra equation. The asymptotio behaviour of the 
ALVE Q oan be oharaoterized as follows: For every antisymmetrie matrix A there exists a unique 
subset of states / C {1,..., S'} whose eoneentrations stay away from zero for all times, that is, 

Xi(t) > Const(A, xo) > 0 for all f > 0 and for every i e I. (6) 

The set I is the set of eondensates. All of the other states with indiees inJ = {1,...,S'} — / 
beeome depleted as f —)■ oo, that is, 

Xi{t) —)■ 0 as f —)■ oo for every i e I. (7) 
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The set of eondensates can be determined algebraically from the antisymmetric matrix A and does 
not depend on the initial conditions xq G As-i = {x G | x* > 0 for all i, 

To show this result, the time-dependent entropy Zi)(c||x)(f) of a condensate vector c = 

(ci,... ,c„) G As-i (Cj > 0 for all i and = 1) relative to the trajectory x(f) is consid¬ 

ered (that is, the Kullback-Leibler divergence of x(f) from c), see equation Q. A condensate 
vector is defined via the properties (see Fig. ??b): 

Cj > 0 and (Ac)* = 0 for f G /, and (8) 

Ci = 0 and (Ac)* < 0 for f G 7 . (9) 

Such a vector can always be found for an antisymmetric matrixpl. Notably, the index set / is unique 
although more than one condensate vector may exist. 


Considering the time derivative of the relative entropy Zi)(c||x)(f) and employing equa¬ 
tions Q and @ yields: 

j ■S' Q S’ s 

—T)(c||x)(f) = - = -^Ci(Ax)i = ^(Ac)iXi = ^(Ac)iXi . (10) 

Cl b Xo _ 

7=1 7=1 7=1 i^I 

Since (Ac) / < 0 and X > 0, it follows that c)ti9(c| |x)(f) < 0. Therefore, the relative entropy 
i9(c| |x) is a Lyapunov function if c is chosen in accordance with equations @ and Q. Moreover, 
i9(c| |x) is bounded from above by i9(c| |x)(0) and from below by zero for all times. This can be 
seen from the definition of D, and from integration of equation ( [T^ (using that (Ac)/ < 0 and 

X > 0): 


0 < i9(c||x)(t) = i9(c||x)(0) -f 



^(Ac)iXi(s) < L>(c||x)(0) . 
iel 


( 11 ) 
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From the definition of the relative entropy in equation Q, it follows that every concentration 
Xi with i G / remains larger than a positive constant, that is, Xi{t) > Const{A,xo) > 0 for all 
times t (if Xi{t) —)■ 0 for i G /, it follows that D ^ oo, which contradicts the boundedness of D). 


Furthermore, equation ([TT]) implies that. 


- ds {Ac)iXi{s) < - ds ^(y4c)ja;i(s) < D(c||x)(0) , 


( 12 ) 


is/ 


for every i E I and for all t. Therefore, concentration Xi is integrable for every i E I {Xi E 
L^(0, oo)) with the bound: 


0 < / ds Xi{s) < = Const{A, xq) for every i e I. 

Jo [Acji 


(13) 


Since the derivative of the concentrations is bounded from above and below, \dtXi\ = \xi{Ax)i\ < 
II (Ax) II oo < IIAII 00^00 < Const (A), one concludes that Xi is uniformly continuous (||A||oo^.oo 
denotes the operator norm of A induced by the maximum norm on M*^). Together with the integra- 
bility ( fT3| ), it follows that states with indices in I become depleted as f —)■ oo, that is, Xi{t) —0 
for i G /. 


In conclusion, given an antisymmetric matrix A = R — RA via a rate constant matrix R = 
one finds a condensate vector c that satisfies inequalities The index set /, for which 

entries of c are positive, represents condensates. The index set I, for which entries of c are zero, 
represents states that become depleted. Moreover, equation ( [T0| ) implies that the relative entropy 
becomes a conserved quantity in the subsystem of condensates (Supplementary Fig. 1). 
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Temporal average of condensate concentrations. The ALVE @ is solved implieitly by, 


x,{t) = , 


with the time average of the trajeetory (x)t defined as: 


1 /■* 

(x)t = - ds x(s) 
t In 


( 14 ) 


(15) 


It is shown above that 0 < Const{A,XQ) < Xi{t) < 1 holds for the states that beeome 
eondensates (i G I). By rearranging equation ( [T4| ), one thus obtains: 


|(4l(x)t)i| < - 


log 


Xijt) 

Xi(0) 


< 


Const{A, xo) 


for alH G / 


(16) 


Note that Const is used to denote arbitrary positive, time-independent eonstants. Therefore, the 


right hand side of equation ( [161 ) vanishes for f —)■ cx). On the other hand, Xi is integrable for i G / 
(equation ([T3|)). Thus, the eorresponding eomponent of the time average eonverges to zero. 


, , Const(A,xo) . - 

{Xi)t < -)■ 0 as f —)■ oo for every i e I . 


(17) 


Henee, the distanee of the time average (x)t to the kernel of the antisymmetrie submatrix A^ 
eonverges to zero (the submatrix A^ eorresponds to the system of eondensates with indiees in J). 


Structure of a generic antisymmetric matrix. For systems with an even number of states S, 
the antisymmetrie matrix A = R — BA generieally has a trivial kernel, whereas for systems with 
an odd number of states, the kernel of A is generieally one-dimensional. A higher dimensional 
kernel of A only oeeurs if the matrix entries are tuned^^lMlZSI. As a eonsequenee, when all of the 
entries above the diagonal of A are, for example, randomly drawn from a eontinuous probability 


18 










distribution (for example from a Gaussian distribution), all 2“^ submatriees of A have a kernel with 
dimension of less than or equal to one. 

The projeetion of x G to the subspaee C for an arbitrary index set J C 
{1,..., 5} is defined as xj := Pjx := {xj)j(zj. In other words, entries of xj are zero for in- 
diees in the eomplement J. In the following, the short notation A'^ := PjAPj is also used (see 
above). Furthermore, the set of antisymmetrie matriees whose submatriees have a kernel with 
dimension less than or equal to one is defined: 

:= {A G I A is antisymmetrie and dimker A"^ < 1 for all J C {1,..., S'}} . (18) 

The eomplement has measure zero with respeet to the flat measure dA on antisymmetrie matriees 
(the translation invariant measure, whieh is sigma-finite and not trivial). 

For antisymmetrie matriees A G fl, the kernel ean be eharaeterized as follow ^^3^. If the 
number of states S is even, the kernel of A is trivial: ker A = {0}. If the number of states is odd, 
the kernel is one-dimensional: ker A = {v}. This kernel element ean be eomputed analytieally in 
terms of Pfaffians of submatriees of A: 

V=(Pf(/li),-Pf(/l2),...,Pf(/l5)) . (19) 

The submatrix Aj^ G denotes the matrix for whieh the k-th eolumn and row are 

removed from A. 

For antisymmetrie matriees A G fl, the normalized eondensate veetor c with G = 1 ^nd 
with properties @ is unique. The latter follows from A^c = 0 (equations @ and ([TS])). There¬ 
fore, the eondensate veetor is the unique kernel veetor of the subsystem of eondensates whose 
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interactions are characterized by the matrix . Furthermore, / contains an odd number of el¬ 
ements. To determine the condensate vector for A G one can proceed as follows. For each 
odd-dimensional submatrix with / C {!,...,5'}, one computes the kernel element v accord¬ 
ing to equation ( [79] ) and defines the vector w G by setting w/ = v and wj = 0. There exists 
exactly one set / for which (Aw)/ < 0. The corresponding vector w is the unique condensate 
vector upon normalization. 


Temporal average of condensate concentrations (generic case). It was shown above that the 
temporal average of condensate concentrations (x)^ converges to a non-negative kernel element of 
the antisymmetric matrix A^. In the generic case, the condensate vector c is the unique kernel ele¬ 
ment of A^ upon normalization. Therefore, positive entries of c represent the asymptotic temporal 
average of condensate concentrations. 


(x)^-c||oo < 


Const{A, xo) 


—)■ 0 as f —)■ cxD 


( 20 ) 


Exponentially fast depletion of states (generic case). Upon inserting equation ( [20| ) into the im¬ 
plicit solution ( [T^ of the ALVE, the exponentially fast depletion of states with i G / can be seen 
as follows (note that (Ac )i < 0 according to the choice of the condensate vector in equations @ 
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and @): 


and analogously, 


Xi{t) = 

(21) 


(22) 


(23) 

= Consf(A, , 

(24) 

Xi{t) = > Const{A,^o)e*<^^^^ . 

(25) 


Therefore, condensate selection occurs exponentially fast at depletion rate |(y4c)i|. The dynamics 
of cases for non-generic antisymmetric matrices are discussed in Supplementary Note 2. 


Linear programming algorithm. For the numeric computation of condensate vectors c, a finite 
threshold 5 > 0 was introduced into the inequalities Ac < 0 and c — Ac > 5 > 0. Its 
value was set to 5 = 1 by rescaling of c. Numerical solution of the inequalities was performed 
by using the IBM ILOG CPLEX Optimization Studio 12.5 and its interface to the C-i-i- language. 
The software Mathematica 9.0 from Wolfram Research was also found to be applicable. Further 
information on the calibration of the linear programming algorithm and a simplified Mathematica 
algorithm are provided in Supplementary Note 3. 
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Antisymmetric Condensate Depietion 

matrix A vector c rates Ac 


Figure 1. Condensation into multiple states due to particle transitions between states and 
mathematics of condensate selection, (a) With respect to condensation in an incoherently driven- 
dissipative quantum system, each bowl represents a state Ei that is occupied by iVj non-interacting 
bosons (filled circles). If indicated by an arrow, bosons may undergo transitions from state Ej to 
state Ei at a rate = rij{Ni + l)Nj, with rate constant Vij. In the language of evolutionary 
game theory, the figure depicts the interaction of iV* agents (filled circles) playing strategies Ei 
(bowls). An agent playing strategy Ej adopts strategy Ei at a rate Ti^^j = rijNiNj. The above rate 
of bosonic condensate selection is recovered if agents may also spontaneously mutate from Ej to 
Ei at a rate r^j. (b) A condensate vector c for an antisymmetric matrix A has two properties: its 
entries are positive for indices for which Ac is zero, and they are zero for indices for which Ac is 
negative (“-” signifies the antisymmetry of matrix A). Temporal evolution of the relative entropy 
of the condensate vector to the state concentrations under the ALVE Q relates positive entries of 
the condensate vector to condensates, and its zero entries to depleted states. Generically, positive 
entries of c represent the asymptotic temporal average of oscillating condensate concentrations 
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according to the ALVE Q, and negative entries of Ac represent depletion rates. 
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Figure 2. Fragmentation of an exemplary system into multiple condensates with oscillat¬ 
ing state concentrations, (a) Randomly sampled network of 50 states. Disks represent states. 
An arrow from state Ej to state Ei represents an effeetive rate eonstant Uij = Vij — r^i = 1 (a 
missing arrow indieates a forbidden transition with aij = 0). Computation of a eondensate vee- 
tor c predieted relaxation into ten isolated eondensates (yellow), one interaeting subsystem with 
six eondensates (blue), and two roek-paper-seissors (RPS) eyeles (red and green). All other states 
beeome depleted. The eomplete network also eomprises RPS eyeles of whieh some states beeome 
depleted. Knowledge of the network topology alone is thus insuffieient to determine eondensates. 
(b) Temporal evolution of state eoneentrations x* (logarithmie seale). Colours in aeeordanee with 
(a). Numerieal integration of the ALVE @ eonfirmed the seleetion of states based on the eonden¬ 
sate veetor c. Subsystems with six (blue) and three (red and green) eondensates exhibit oseillations 
of eoneentrations with non-vanishing partiele flow. Depletion of states oeeurs exponentially fast. 
Identifying eondensates from a eondensate veetor c is more reliable than through numerieal inte¬ 
gration: The eoneentration of the state assoeiated to the purple disk in (a) deeays exponentially to 
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a concentration of 1.5 ■ 10 ^ before recovering transiently. Numerical integration cannot rule out 
permanent recovery at later times. Supplementary Fig. 2 demonstrates such a case. 
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Figure 3. Dependence of the number of condensates on the connectivity of states in random 
networks, (a-f) Measured probability p of finding a particular number of condensates for a system 
with S = 100 states and with connectivity C (5 ■ 10® systems analysed per histogram). The 
connectivity specifies the percentage of states between which transitions of particles occur with a 
non-zero effective rate constant Ujj = rij — rji. Effective rate constants Uij were sampled from 
a Gaussian distribution (zero mean, unit variance), (a) At full connectivity, the distribution is 
pseudo-binomial with only odd numbers of condensates (G = 1; light blue bars (b) As the 
connectivity is reduced, even numbers of condensates become possible when systems decouple into 
even numbers of subsystems {C = 0.129; dark blue bars), (c) The distribution exhibits bimodality 
(C = 0.075) and (d) approaches a minimal average number of 40.2 condensates (C = 0.055). (e) 
This average subsequently increases (C = 0.018) because isolated states are trivially selected as 


35 

































condensates (C = 0) as shown in (f). (g) Average number of eondensates per number of states 
(eolour-eoded) plotted against the number of states S and the eonneetivity C (log-log graph in 
inset; > 10"^ systems per data point, see Supplementary Fig. 4 for the reliability of the linear 
programming algorithm). White eireles eorrespond to distributions shown in (a-e). The minimal 
relative number of eondensates eonforms to the power law C ~ 1/S"’' with 7 = 0.998 ± 0.008 
(s.e.m.) and ean be related to the eritieality of random network^. 
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RPS condition 


Attractivity condition 


Figure 4. Conditions for the emergence of a rock-paper-scissors cycle of condensates. Three 
particular states Ei,E 2 , and E 3 (blue, red, and yellow disks) of a network condense into a rock- 
paper-scissors (RPS) cycle if, and only if, two conditions are fulfilled: First, the “RPS condition” 
requires that the rate constants rij between the three states form a RPS network: rj_i ,+1 > 

(indices are counted modulo 3, for example, r ^2 = t^i 2 (framed arrows denote rate constants that 
are larger than rate constants for the respective reverse direction). Differences between these rate 
constants define the entries q = of an admissible condensate vector c. Second, 

the “attractivity condition” requires that the weighted sum of rates from any exterior state E^ 
(purple disks) into the RPS cycle, (framed arrows), is larger than the weighted sum of 

outbound rates, ^j'^kj (black arrows). In other words, the inflow of particles into the RPS 
cycle from any exterior state needs to be greater than the outflow to that state. 
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Field of research 

Entity 

State 

Process 

Dynamics 

Quantum physics^'’ 

Boson 

Quantum state 

Incoherent transition 

Condensation/depletion 

Evolutionary game theory^’’ ™ 

Agent 

Strategy 

Game (+mutation) 

Win/loss 

Population dynamics®^'®^ '*^ 

Individual 

Species 

Competition (+mutation) 

Survival/extinction 

Chemical kinetics^® ®® 

Molecule 

Chemical species 

Reaction (+conversion) 

Production/consumption 

Plasma physics®'*' ®® 

Plasmon 

Jet 

Scattering 

Increase/decrease 


Table 1 : Condensation processes described by the ALVE in different fields of research. The 
ALVE @ governs condensation processes in diverse fields of research. For example, for 
incoherently driven-dissipative bosonic systems, the ALVE describes condensation and 
depletion of states by incoherent transitions of non-interacting bosons. In EOT, the ALVE 
occurs in the context of winning and losing strategies played by agents. 
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